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B DECAYS AND THE HEAVY-QUARK EXPANSION 

MATTHIAS NEUBERT 
Theory Division, CERN, CH-1211 Geneva 23, Switzerland 

We review the theory and phenomenology of heavy-quark symmetry, exclusive 
weak decays of B mesons, inclusive decay rates and lifetimes of b hadrons. 

1 Introduction 

The rich phenomenology of weak decays has always been a source of informa- 
tion about the nature of elementary particle interactions. A long time ago, 
/3- and /i-decay experiments revealed the structure of the effective flavour- 
changing interactions at low momentum transfer. Today, weak decays of 
hadrons containing heavy quarks are employed for tests of the Standard Model 
and measurements of its parameters. In particular, they offer the most direct 
way to determine the weak mixing angles, to test the unitarity of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix, and to explore the physics of CP violation. 
On the other hand, hadronic weak decays also serve as a probe of that part of 
strong-interaction phenomenology which is least understood: the confinement 
of quarks and gluons inside hadrons. 

The structure of weak interactions in the Standard Model is rather simple. 
Flavour-changing decays are mediated by the coupling of the charged current 
Jqq to the M^-boson field: 

rcc^-^J^cW^P^+h.c, (1) 

where 

JcC^ ('^e,'^p.,t^T)l^ A^L + (uL,CL,tL)7^^CKM Sl (2) 

contains the left-handed lepton and quark fields, and 

(Vud Vus Vub \ 
Vcd Vcs Vch (3) 

Vtd Vts VtbJ 

is the CKM matrix. At low energies, the charged-current interaction gives rise 
to local four-fermion couplings of the form 

£eff = -2V2Gi.Jgc4c.^> (4) 



where 



Gf- 



4V2M2 



1.16639(2) GeV" 



(5) 



w 



is the Fermi constant. 
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Figure 1: Examples of leptonic {B — > r Ut), semileptonic (B" — » D+e Ue), and non- 
leptonic {B'^ — » D+7r~) decays of B mesons. 



According to the structure of the charged-current interaction, weak decays 
of hadrons can be divided into three classes: leptonic decays, in which the 
quarks of the decaying hadron annihilate each other and only leptons appear in 
the final state; semileptonic decays, in which both leptons and hadrons appear 
in the final state; and non-leptonic decays, in which the final state consists 
of hadrons only. Representative examples of these three types of decays are 
shown in Fig. n^. The simple quark-line graphs shown in this figure are a 
gross oversimplification, however. In the real world, quarks are confined inside 
hadrons, bound by the exchange of soft gluons. The simplicity of the weak 
interactions is overshadowed by the complexity of the strong interactions. A 
complicated interplay between the weak and strong forces characterizes the 



phenomenology of hadronic weak decays. As an example, a more realistic 
picture of a non-leptonic decay is shown in Fig. 0. 
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Figure 2: More realistic representation of a non-leptonic decay. 



The complexity of strong-interaction effects increases with the number of 
quarks appearing in the final state. Bound-state effects in leptonic decays 
can be lumped into a single parameter (a "decay constant"), while those in 
semileptonic decays are described by invariant form factors, depending on the 
momentum transfer q^ between the hadrons. Approximate symmetries of the 
strong interactions help us to constrain the properties of these form factors. 
For non-leptonic decays, on the other hand, we are still far from having a 
quantitative understanding of strong-interaction effects even in the simplest 
decay modes. 

Over the last decade, a lot of information on heavy-quark decays has been 
collected in experiments at e+e^ storage rings operating at the T(4s) res- 
onance, and more recently at high-energy e^e^ and hadron colliders. This 
has led to a rather detailed knowledge of the flavour sector of the Standard 
Model and many of the parameters associated with it_jThere have been sevieral 
great discoveries in this field, such as B^-B^ mbcingEJa, b ^ u transitions 013, 
and rare decays induced by penguin operatorsu. Yet there is much more to 
come. In particular, the S-factorics at SLAC, KEK, HERA-B and LHC-B will 
provide a wealth of new results within the coming years. 

The experimental progress in heavy-fiavour physics has been accompa- 
nied by a significant progress in theory, which was related to the discovery 
of heavy-quark symmetry, the development of the heavy-quark effective the- 
ory, and the establishment of the heavy-quark expansion. The excitement 



about these developments rests upon the fact that they aUow (some) model- 
independent predictions in an area in which "progress" in theory often meant 
nothing more than the construction of a new model, which could be used to 
estimate some strong-interaction hadronic matrix elements. In section g, we 
explain the physical picture behind heavy-quark symmetry and discuss the con- 
struction, as well as simple applications, of the heavy-quark effective theory. 
Section || deals with applications of these concepts to exclusive weak decays 
of B mesons. Applications of the heavy-quark expansion to the description of 
inclusive decay rates and lifetimes of b hadrons are the topics of section |4[ 

2 Heavy-Quark Symmetry 

This paction provides an introduction to thCpideas of heavy-quark symme- 
try I3"ll3 and the heavy-quark effective theory Il3~c^, which provide the mod- 
ern theoretical framework for the description of the properties and decays of 
hadrons containing a heavy quark. For a more detailed description of this 
subject, the reader is referred to the review articles in Refs. 24-30. 

2.1 The Physical Picture 

There are several reasons why the strong interactions of hadrons containing 
heavy quarks are easier to understand than those of hadrons containing only 
light quarks. The first is asymptotic freedom, the fact that the effective cou- 
pling constant of QCD becomes weak in processes with a large-momentum 
transfer, corresponding to interactions at short distance scales EilS. At large 
distances, on the other hand, the coupling becomes strong, leading to non- 
perturbative phenomena such as the confinement of quarks and gluons on a 
length scale i?had ~ 1/Aqcd ^ 1 fm, which determines the size of hadrons. 
Roughly speaking, Aqcd ^ 0.2 GeV is the energy scale that separates the 
regions of large and small coupling constant. When the mass of a quark Q 
is much larger than this scale, mq 3> Aqcd, it is called a heavy quark. The 
quarks of the Standard Model fall naturally into two classes: up, down and 
strange are light quarks, whereas charm, bottom and top are heavy quarks q 
For heavy quarks, the effective coupling constant as{mQ) is small, implying 
that on length scales comparable to the Compton wavelength Xq ^ ^/rriQ 
the strong interactions are perturbative and much like the electromagnetic 
interactions. In fact, the quarkonium systems (QQ), whose size is of order 
Xg/asimQ) <C -Rhad, are very much hydrogen- like. 



"Ironically, the top quark is of no relevance to our discussion here, since it is too heavy 
to form hadronic bound states before it decays. 



Systems composed of a heavy quark and other hght constituents are more 
comphcated. The size of such systems is determined by i?had, and the typi- 
cal momenta exchanged between the heavy and Ught constituents are of order 
Aqcd- The heavy quark is surrounded by a most comphcated, strongly inter- 
acting cloud of light quarks, antiquarks and gluons. In this case it is the fact 
that Ag <C -Rhad, i-e. that the Compton wavelength of the heavy quark is much 
smaller than the size of the hadron, which leads to simplifications. To resolve 
the quantum numbers of the heavy quark would require a hard probe; the soft 
gluons exchanged between the heavy quark and the light constituents can only 
resolve distances much larger than Ag. Therefore, the light degrees of freedom 
are blind to the flavour (mass) and spin orientation of the heavy quark. They 
experience only its colour field, which extends over large distances because of 
confinement. In the rest frame of the heavy quark, it is in fact only the elec- 
tric colour field that is important; relativistic effects such as colour magnetism 
vanish as mg — *■ cx). Since the heavy-quark spin participates in interactions 
only through such relativistic effects, it decouples. That the heavy-quark mass 
becomes irrelevant can be seen as follows: As ruQ — > oo, the heavy quark and 
the hadron that contains it have the same velocity. In the rest frame of the 
hadron, the heavy quark is at rest, too. The wave function of the light con- 
stituents follows from a solution of the field equations of QCD subject to the 
boundary condition of a static triplet source of colour at the location of the 
heavy quark. This boundary condition is independent of mg, and so is the 
solution for the configuration of the light constituents. 

It follows that, in the limit mg — > oo, hadronic systems which differ only 
in the fiavour or spin quantum numbers of J;h£.heavy quark have the same con- 
figuration of their light degrees of freedomDHiS. Although this observation still 
does not allow us to calculate what this configuration is, it provides relations 
between the properties of such particles as the heavy mesons B, D, B* and D*, 
or the heavy baryons Ab and Ac (to the extent that corrections to the infinite 
quark-mass limit are small in these systems). These relations result from some 
approximate symmetries of the effective strong interactions of heavy quarks 
at low energies. The configuration of light degrees of freedom in a hadron 
containing a single heavy quark with velocity v does not change if this quark 
is replaced by another heavy quark with different flavour or spin, but with the 
same velocity. Both heavy quarks lead to the same static colour field. For N/^ 
heavy-quark fiavours, there is thus an SU(2iV/i) spin-flavour symmetry group, 
under which the effective strong interactions are invariant. These symmetries 
are in close correspondence to familiar properties of atoms. The fiavour sym- 
metry is analogous to the fact that different isotopes have the same chemistry, 
since to good approximation the wave function of the electrons is independent 



of the mass of the nucleus. The electrons only see the total nuclear charge. 
The spin symmetry is analogous to the fact that the hyperfine levels in atoms 
are nearly degenerate. The nuclear spin decouples in the limit me/mM —^ 0. 

Heavy-quark symmetry is an approximate symmetry, and corrections arise 
since the quark masses are not infinite. In many respects, it is complementary 
to chiral symmetry, which arises in the opposite limit of small quark masses. 
There is an important distinction, however. Whereas chiral symmetry is a 
symmetry of the QCD Lagrangian in the limit of vanishing quark masses, 
heavy-quark symmetry is not a symmetry of the Lagrangian (not even an ap- 
proximate one), but rather a symmetry of an effective theory, which is a good 
approximation of QCD in a certain kinematic region. It is realized only in sys- 
tems in which a heavy quark interacts predominantly by the exchange of soft 
gluons. In such systems the heavy quark is almost on-shell; its momentum fluc- 
tuates around the mass shell by an amount of order Aqcd- The corresponding 
fluctuations in the velocity of the heavy quark vanish as Aqcd /niQ —^ 0. The 
velocity Jaecomes a conserved quantity and is no longer a dynamical degree of 
freedoml23. Nevertheless, results derived on the basis of heavy-quark symme- 
try are model-independent consequences of QCD in a well-defined limit. The 
symmetry-breaking corrections can be studied in a systematic way. To this 
end, it is however necessary to cast the QCD Lagrangian for a heavy quark, 

CQ=Q{tp-mQ)Q, (6) 

into a form suitable for taking the limit mq — > oo. 

2.2 Heavy- Quark Effective Theory 

The effects of a very heavy particle often become irrelevant at low energies. It 
is then useful to construct a low-energy effective theory, in which this heavy 
particle no longer appears. Eventually, this effective theory will be easier to 
deal with than the full theory. A familiar example is Fermi's theory of the 
weak interactions. For the description of the weak decays of hadrons, the 
weak interactions can be approximated by point-like four-fermion couplings, 
governed by a dimensionful coupling constant Gp [cf . (|j)] . The effects of the 
intermediate vector bosons, W and Z , can only be resolved at energies much 
larger than the hadron masses. 

The process of pemoying the degrees of freedom of a heavy particle involves 
the following steps Lj^Ej: one first identifies the heavy-particle fields and "in- 
tegrates them out" in the generating functional of the Green functions of the 
theory. This is possible since at low energies the heavy particle does not appear 
as an external state. However, although the action of the full theory is usually 



a local one, what results after this first step is a non-local effective action. The 
non-locality is related to the fact that in the full theory the heavy particle 
with mass M can appear in virtual processes and propagate over a short but 
finite distance Aa; ^ 1/M. Thus, a second step is required to obtain a local 
effective Lagrangian: the non-local effective action is rewritteia_as_an infinite 
series of local terms in an Operator Product Expansion (OPE) EjO. Roughly 
speaking, this corresponds to an expansion in powers of 1/M. It is in this step 
that the short- and long-distance physics is disentangled. The long-distance 
physics corresponds to interactions at low energies and is the same in the full 
and the effective theory. But short-distance effects arising from quantum cor- 
rections involving large virtual momenta (of order M) are not reproduced in 
the effective theory, once the heavy particle has been integrated out. In a third 
step, they have to be added in a perturbative way using renormalization-group 
techniques. These short-distance effects lead to a renormalization of the coef- 
ficients of the local operators in the effective Lagrangian. An example is the 
effective Lagrangian for non-leptonic weak decays, in which radiative correc- 
tions from hard gluons with virtual momenta in the range between myy and 
some renormalization scale fj, ^ 1 GeV give rise^tp- Wilson coefficients, which 
renormalize the local four-fermion interactions l3lj. 

The heavy-quark effective theory (HQET) is constructed to provide a sim- 
plified description of processes where a heavy quark interaitSpSSiith light degrees 
of freedom predominantly by the exchange of soft gluons tjc3. Clearly, mg is 
the high-energy scale in this case, and Aqcd is the scale of the hadronic physics 
we are interested in. The situation is illustrated in Fig. 0. At short distances, 
i.e. for energy scales larger than the heavy-quark mass, the physics is pertur- 
bative and described by conventional QCD. For mass scales much below the 
heavy-quark mass, the physics is complicated and non-perturbative because 
of confinement. Our goal is to obtain a simplified description in this region 
using an effective field theory. To separate short- and long-distance effects, we 
introduce a separation scale fj, such that Aqcd ^ /i <C mq. The HQET will 
be constructed in such a way that it is identical to QCD in the long-distance 
region, i.e. for scales below /i. In the short- distance region, the effective the- 
ory is incomplete, since some high-momentum modes have been integrated out 
from the full theory. The fact that the physics must be independent of the 
arbitrary scale fj, allows us to derive renormalization-group equations, which 
can be employed to deal with the short-distance effects in an efficient way. 

Compared with most effective theories, in which the degrees of freedom 
of a heavy particle are removed completely from the low-energy theory, the 
HQET is special in that its purpose is to describe the properties and decays of 
hadrons which do contain a heavy quark. Hence, it is not possible to remove 
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Figure 3: Philosophy of the heavy-quark eflfeetive theory. 



the heavy quark completely from the effective theory. What is possible is to 
integrate out the "small components" in the full heavy-quark spinor, which 
describe the fluctuations around the mass shell. 

The starting point in the construction of the HQET is the observation that 
a heavy quark bound inside a hadron moves more or less with the hadron's 
velocity v, and is almost on-shell. Its momentum can be written as 



Pq 



mQV^ 



k" . 



(7) 



where the components of the so-called residual momentum k arc much smaller 
than rUQ. Note that w is a four- velocity, so that w^ = 1. Interactions of the 
heavy quark with light degrees of freedom change the residual momentum by 
an amount of order Afc ^ Aqcd , but the corresponding changes in the heavy- 
quark velocity vanish as Aqcd/tiq -^ 0. In this situation, it is appropriate to 
introduce large- and small-component fields, h^ and H^, by 

h, {x) = e™^'^-" P+ Q{x) , H, {x) - e^"'?'^-" P_ Q{x) , (8) 

where P-f and P_ are projection operators defined as 



It follows that 

Q{x) = e-^'"«"-=^ [K{x) + H,{x)] . (10) 

Because of the projection operators, the new fields satisfy fhy = hy and f Hy = 
—Hy. In the rest frame, i.e. for u^ = (1,0,0,0), hy corresponds to the upper 
two components of Q, while Hy corresponds to the lower ones. Whereas hy 
annihilates a heavy quark with velocity v, Hy creates a heavy antiquark with 
velocity v. 

In terms of the new fields, the QCD Lagrangian (g) for a heavy quark 
takes the form 

^Q ~ hy iv-Dhy — Hy iiv-D + 2toq) Hy + hy ip±Hy + Hy ipj_hy , (11) 

where D'^ = _D^ — v^^v-Dis orthogonal to the heavy-quark velocity: v ■ D± — 0. 
In the rest frame, Dj = (0, D ) contains the spatial components of the covariant 
derivative. From (|l]), it is apparent that hy describes massless degrees of 
freedom, whereas Hy corresponds to fluctuations with twice the heavy-quark 
mass. These are the heavy degrees of freedom that will be eliminated in the 
construction of the effective theory. The fields are mixed by the presence of the 
third and fourth terms, which describe pair creation or annihilation of heavy 
quarks and antiquarks. As shown in the first diagram in Fig. ^ in a virtual 
process, a heavy quark propagating forward in time can turn into an antiquark 
propagating backward in time, and then turn back into a quark. The energy of 
the intermediate quantum state hhH is larger than the energy of the incoming 
heavy quark by at least 2mQ. Because of this large energy gap, the virtual 
quantum fluctuation can only propagate over a short distance Ax ~ l/rriQ. 
On hadronic scales set by i?had = 1/Aqcd, the process essentially looks like a 
local interaction of the form 

hyip±- ip±hy, (12) 

2mQ 

where we have simply replaced the propagator for Hy by l/2mQ. A more 
correct treatment is to integrate out the small-component field Hy, thereby 
deriving a non-local effective action for the large-component field hy, which 
can then be expanded in terms of local operators. Before doing this, let us 
mention a second type of virtual corrections involving pair creation, namely 
heavy-quark loops. An example is shown in the second diagram in Fig. m 
Heavy-quark loops cannot be described in terms of the effective fields hy and 
Hy, since the quark velocities inside a loop are not conserved and are in no 
way related to hadron velocities. However, such short-distance processes are 
proportional to the small coupling constant asimq) and can be calculated in 

9 



perturbation theory. They lead to corrections that are added onto the low- 
energy effective theory in the renormalization procedure. 





Figure 4: Virtual fluctuations involving pair creation of heavy quarks. Time flows to the 

right. 

On a classical level, the heavy degrees of freedom represented by Hy can 
be eliminated using the equation of motion. Taking the variation of the La- 
grangian with respect to the field H^, we obtain 

{ivD + 2mQ)H^^ipj_K. (13) 

This equation can formally be solved to give 

Hv^l^ ^-—f,^P±K, (14) 

zniQ + iv-D 

showing that the small-component field H^ is indeed of order 1/toq. We can 
now insert this solution into (^l]) to obtain the "non-local effective Lagrangian" 

Ccff = Kiv-Dhy + hyip±- ; — —ip±hy. (15) 

ZniQ + tv-D 

Clearly, the second term corresponds to the first class of virtual processes 
shown in Fig. ^. 

It is possible to derive this Lagrangian in a more elegant way by manipulat- 
ing thajgenerating functional for QCD Green functions containing heavy-quark 
fieldslH3. To this end, one starts from the field redefinition (nfl) and couples 
the large-component fields h^ to external sources p^. Green functions with an 
arbitrary number of hy fields can be constructed by taking derivatives with 
respect to py . No sources are needed for the heavy degrees of freedom repre- 
sented by Hy. The functional integral over these fields is Gaussian and can be 
performed explicitly, leading to the effective action 

'S'off = / d^'x £cfr - « In A , (16) 

with CcS as given in (nq) . The appearance of the logarithm of the determinant 

A = exp ( - Tr In [2mQ + ivD- irj] j (17) 

10 



is a quantum effect not present in the classical derivation presented above. 
However, in this case the determinant can be regulated in a gauge-invariant 
way, and by choosingJjhe gauge v ■ A = one can show that In A is just an 
irrelevant constant ac^. 

Because of the phase factor in (|l^) , the x dependence of the effective heavy- 
quark field hy is weak. In momentum space, derivatives acting on hy produce 
powers of the residual momentum k, which is much smaller than toq. Hence, 
the non-local effective Lagrangian (051) allows for a derivative expansion: 



CcS — hyiv-Dhy 



_. oo 

y^ hyil 



n=0 



iv ■ D 



2m, 



Q 



ip±hy 



(18) 



Taking into account that hy contains a P+ projection operator, and using the 
identity 



P+ip^tp^P+ = P+ 



i.D^Y + ^a^.G^^'' 



P+ 



where i[D^, D^] = gs G^"^ is the gluon field-strength tensor, one finds that 



(19) 
0i 



1 



gs 



C^ff = hy IVD hy + -^- hy {l D ^f hy + -f^ hy fj^^ C" hy + O {l / 171%) . (20) 

In the limit mg —^ cxd, only the first terms remains: 

Coo = hyiv-Dhy . (21) 

This is the effective Lagrangian of the HQET. It gives rise to the Feynman 
rules shown in Fig. |^. 
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Figure 5: Feynman rules of the HQET (i,j and a are colour indices). A heavy quark with 
velocity v is represented by a double line. The residual momentum k is defined in (m). 

Let us take a moment to study the symmetries of this LagrangianEH Since 
there appear no Dirac matrices, interactions of the heavy quark with gluons 
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leave its spin unchanged. Associated with this is an SU(2) symmetry group, 
under which Loo is invariant. The action of this symmetry on the heavy-quark 
fields becomes most transparent in the rest frame, where the generators S"* of 
SU(2) can be chosen as 



^'=2(^0 ;.j; \S\S^\^^,^^-S^. (22) 

Here cr* are the Pauli matrices. An infinitesimal SU(2) transformation Kt, — > 
(1 + it- S) hy leaves the Lagrangian invariant: 

SCoo = hy[iv-D,ie- S]hy — . (23) 

Another symmetry of the HQET arises since the mass of the heavy quark does 
not appear in the effective Lagrangian. For N^ heavy quarks moving at the 
same velocity, eq. (|2l| ) can be extended by writing 

Nh _ 

£o. = J2~^liv-Dhl. (24) 

1=1 

This is invariant under rotations in flavour space. When combined with the 
spin symmetry, the symmetry groTm.is, promoted to SV{2Nh). This is the 
heavy-quark spin-flavour symmetry E3E3. Its physical content is that, in the 
limit TTiQ — > 00, the strong interactions of a heavy quark become independent 
of its mass and spin. 

Consider now the operators appearing at order l/mq in the effective La- 
grangian (|20|). They are easiest to identify in the rest frame. The first operator, 

Okin = 7^ K {iD^)^ hy -^ --^ hy {iD f hy , (25) 

2toq 2mQ 

is the gauge-covariant extension of the kinetic energy arising from the residual 
motion of the heavy quark. The second operator is the non-Abelian analogue 
of the Pauli interaction, which describes the colour-magnetic coupling of the 
heavy-quark spin to the gluon field: 

Cmag — -j hy au,i, G^" hy — > hy S-Bchy . (26) 

AniQ TTiQ 

Here S is the spin operator defined in (|2^), and B* — —^e^^^G^'^ are the 
components of the colour-magnetic field. The chromo-magnetic interaction is 
a relativistic effect, which scales like l/mq. This is the origin of the heavy- 
quark spin symmetry. 

12 



2.3 The Residual Mass Term and the Definition of the Heavy-Quark Mass 

The choice of the expansion parameter in the HQET, i.e. the definition of the 
heavy-quark mass mg, deserves some comments. In the derivation presented 
earher in this section, we chose mq to be the "mass in the Lagrangian" , and 
using this parameter in the phase redefinition in (M) we obtained the effective 
Lagrangian (|2l]), in which the heavy-quark mass no longer appears. However, 
this treatment has its subtleties. The symmetries of the HQET allow a "resid- 
ual mass" dm for the heavy quark, provided that dm is of order Aqqd and is 
the same for all heavy-quark flavours. Even if we arrange that such a mass 
term is not present at the tree level, it will in general be induced by quantum 
corrections. (This is unavoidable if the theory is regulated with a dimensionful 
cutoff.) Therefore, instead of (|2l| ) we should write the effective Lagrangian in 
the more general formCj 

'Coo — hyiv-D hy — Sm hyhy . (27) 

If we redefine the expansion parameter according to mq -^ mq + A?77,, the 
residual mass changes in the opposite way: 5m — > 5m — Arn. This implies that 
there is a unique choice of the expansion parameter mq such that 5m = 0. 
Requiring 5m = 0, as it is usually done implicitly in the HQET, defineSpa,heavy- 
quark mass, which in perturbation theory coincides with the pole massEJ. This, 
in turn, defines for each heavy hadron Hq a parameter A (sometimes called 
the "binding energy" ) through 

A = [mHQ - mq) . (28) 

mQ — >oo 

If one prefers to work with another choice of the expansion parameter, the 
values of non-pcrturbative parameters such as A change, but at the same time 
one has to include the residual mass term in the HQET Lagrangian. It can 
be shown that the various parameters depending on the definition of mq en- 
ter the predictions for physical quantities in^iiich a way that the results are 
independent of the particular choice adopted l3. 

There is one more subtlety hidden in the above discussion. The quantities 
TTiQ, A and 5m are non-perturbative parameters of the HQET, which have a 
similar status as the vacuum condensates in QCD phenomenology l3. These 
parameters cannot be defined unambiguously in perturbation theory. The rea- 
son lies in the divergent behaviour of perturbative expansions in large orders, 
which is associated with the existence oLaioftularities along the real axis in the 
Borel plane, the so-called renormalons ea^Ej. For instance, the perturbation 
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series which relates the pole mass toq of a heavy quark to its bare mass, 
"^Q = '"q '° (l + ci asimq) + C2 a^lmg) + . . . + c„ <(mQ) + . . . | , (29) 

contains numerical coefficients c„ that grow, as n\ for large n, rendering the 
series divergent and not Borel summable aCJ. The best one can achieve is 
to truncate the perturbation series at the minimal term, but this leads to an 
unavoidable arbitrariness of order Attiq ~ Aqcd (the size of the minimal term) 
in the value of the pole mass. This observation, which at first sight seems a 
serious problem for QCD phenomenology, should not come as a surprise. We 
know that because of confinement quarks do not appear as physical states 
in nature. Hence, there is no unique way to define their on-shell properties 
such as a pole mass. In view of this, it is actually remarkable that QCD 
perturbation theory "knows" about its incompleteness and indicates, through 
the appearance of renormalon singularities, the presence of non-perturbative 
effects. One must first specify a scheme how to truncate the QCD perturbation 
series before non-perturbative statements such as 6m — become meaningful, 
and hence before non-perturbative parameters such as mg and A become well- 
defined quantities. The actual values of these parameters will depend on this 
scheme. 

We stress that the "renormalon ambiguities" are not a conceptual problem 
for the heavy-quark expansion. In fact, it can be shown quite generaUv that 
these ambiguities cancel in all predictions for physical observablesE3"L3. The 
way the cancellations occur is intricate, however. The generic structure of the 
heavy-quark expansion for an observable is of the form: 

Observable - Clasimn)] f 1 + + ...), (30) 

V rriQ J 

where C[as{mQ)] represents a perturbative coefficient function, and A is a 
dimensionful non-perturbative parameter. The truncation of the perturba- 
tion series defining the coefficient function leads to an arbitrariness of order 
Aqcd/^t-q, which cancels against a corresponding arbitrariness of order Aqcd 
in the definition of the non-perturbative parameter A. 

The renormalon problem poses itself when one imagines to apply pertur- 
bation theory in very high orders. In practise, the perturbative coefficients are 
known to finite order in a^ (at best to two-loop accuracy), and to be consistent 
one should use them in connection with the pole mass (and A etc.) defined to 
the same order. 
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2-4 Spectroscopic Implications 

The spin-flavour symmetry leads to many interesting relations between the 
properties of hadrons containing a heavy niiark. The most direct consequences 
concern the spectroscopy of such states lJ. In the limit tuq —> oo, the spin 
of the heavy quark and the total angular momentum j of the light degrees 
of freedom are separately conserved by the strong interactions. Because of 
heavy-quark symmetry, the dynamics is independent of the spin and mass 
of the heavy quark. Hadronic states can thus be classified by the qjuantum 
numbers (flavour, spin, parity, etc.) of their light degrees of freedom eD. The 
spin symmetry predicts that, for fixed j ^ 0, there is a doublet of degenerate 
states with total spin J — j zt ^. The flavour symmetry relates the properties 
of states with different heavy-quark flavour. 

In general, the mass of a hadron Hq containing a heavy quark Q obeys an 
expansion of the form 



mHQ^mQ + A+^— + 0{l/ml). (31) 

The parameter A represents contributions arisingJ'roin terms in the Lagrangian 
that are independent of the heavy-quark mass l3, whereas the quantity Arn^ 
originates from the terms of order 1/toq in the effective Lagrangian of the 
HQET. For the ground-state pseudoscalar and vector mesons, one can parame- 
trize the contributions from the kinetic energy and the chromo-njLa|gnetic in- 
teraction in terms of two quantities Ai and A2, in such a way thatEJ 



Am^ = -Ai + 2 



J(J+l)-f A2. (32) 



The hadronic parameters A, Ai and A2 are independent of niQ. They charac- 
terize the properties of the light constituents. 

Consider, as a first example, the SU(3) mass splittings for heavy mesons. 
The heavy-quark expansion predicts that 

TUBS - ™Bd == ^s-^d + 0{l/mt) , 

rriDs - "iDd = As - Ad + 0(l/mc) , (33) 

where we have indicated that the value of the parameter A depends on the 
flavour of the light quark. Thus, to the extent that the charm and bottom 
quarks can both be considered sufficiently heavy, the mass splittings should 
be similar in the two systems. This prediction is confirmed experimentally, 
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rriBs - niBa = (90 ± 3) MeV , 

mos - moa = (99 ± 1) MeV . (34) 

As a second example, consider the spin splittings between the ground-state 
pseudoscalar (J — 0) and vector (J = 1) mesons, which are the members of 
the spin-doublet with j = ^- From (|3|) and (P2|), it follows that 

m%, - m% = 4A2 + 0{l/mb) , 

ml, - m,l = 4A2 + 0(l/m,) . (35) 

The data are compatible with this: 

m%. - m| « 0.49 GeV^ , 

ml, -ml K 0.55 GeV^ . (36) 

Assuming that the B system is close to the heavy-quark limit, we obtain the 
value 

A2 ~ 0.12 GeV^ (37) 

for one of the hadronic parameters in (|32| ) . This quantity plays an important 
role in the phenomenology of inclusive decays of heavy hadrons. 

A third example is provided by the mass splittings between the ground- 
state mesons and baryons containing a heavy quark. The HQET predicts that 

"lAi, - mB = Abaryon ^ A,„cson + 0{l/mb) , 

TTiA, - mo = Abaryon " A,„oson + 0{l/mc) . (38) 

This is again consistent with the experimental results 

TOA, - TOB = (346 ± 6) MeV , 

TOA, -mo = (416 ± 1) MeV , (39) 

although in this case the data indicate sizeable symmetry-breaking corrections. 
For the mass of the A;, baryon, we have used the value 

mA, == (5625 ± 6) MeV , (40) 

which is obtained by averaging the resultEd m,\^ = (5639 ± 15) MeV with the 
value mA, = (5623 ± 5 ± 4) MeV reported by the CDF CoUaborationN. The 
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dominant correction to the relations (Bq) comes from the contribution of the 
chromo-magnetic interaction to the masses of the heavy mesons fl which adds 
a term 3A2/2mQ on the right-hand side. Including this term, we obtain the 
refined prediction that the two quantities 

ruA, - niB = (312 ± 6) MeV . 

2mB 

TOA^ -mo- -^-^ = (320 ± 1) MeV (41) 

2mD 

should be close to each other. This is clearly satisfied by the data. 

The mass formula ( |3l| ) can also be used to derive information on the 
heavy-quark masses from the observed hadron masses. Introducing the "spin- 
averaged" meson masses ttTb = \{mB + Stub*) ~ 5.31 GeV and toc = 
J {mo + SniD-) « 1.97 GeV, we find that 

mb-mc = {rriB - mo) \ 1 - ^_ L + 0{\/m%) \ . (42) 

Using theoretical estimates for the parameter Ai, which lie in the range EJ~0 

Ai =-(0.3±0.2) GeV^ (43) 

this relation leads to 

rrib - nic = (3.39 ± 0.03 ± 0.03) GeV , (44) 

where the first error reflects the uncertainty in the value of Ai, and the second 
one takes into account unknown higher-order corrections. The fact that the 
difference {mb — rric) is determined rather precisely becomes important in the 
analysis of inclusive decays of heavy hadrons. 

3 Exclusive Semileptonic Decays 

Semileptonic decays of B mesons have received a lot of attention in recent 
years. The decay channel B -^ D*tv has the largest branching fraction of all 
_B-meson decay modes. From a theoretical point of view, semileptonic decays 
are simple enough to allow for a reliable, quantitative description. The analysis 
of these decays provides much information about the strong forces that bind the 



''Because of the spin symmetry, there is no such contribution to the masses of the Ag 
baryons. 
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quarks and gluons into hadrons. Schematically, a semileptonic decay process 
is shown in Fig. o. The strength of the b ^ c transition vertex is governed 
by the element Vet of the CKM matrix. The parameters of this matrix are 
fundamental parameters of the Standard Model. A primary goal of the study 
of semileptonic decays of B mesons is to extract with high precision the values 
of \Vcb\ and \Vub\- We will now discuss the theoretical basis of such analyses. 



^cbX c 




D. D^. D^ 



Figure 6: Semileptonic decays of B mesons. 



3.1 Weak Decay Form Factors 

Heavy-quark symmetry implies relations between the weak decay form factors 
of heavy mesons, which are_Df particular interest. These relations have been 
derivecLby Isgur and Wise 113, generalizing .ideas developed by Nussinov and 
Wetzel I, and by Voloshin and ShifmanElEl. 

Consider the elastic scattering of a -B meson, B{v) -^ B{v'), induced by a 
vector current coupled to the b quark. Before the action of the current, the light 
degrees of freedom inside the B meson orbit around the heavy quark, which 
acts as a static source of colour. On average, the b quark and the B meson have 
the same velocity v. The action of the current is to replace instantaneously 
(at time t — to) the colour source by one moving at a velocity u', as indicated 
in Fig. 0. li V = v' , nothing happens; the light degrees of freedom do not 
realize that there was a current acting on the heavy quark. If the velocities are 
different, however, the light constituents suddenly find themselves interacting 
with a moving colour source. Soft gluons have to be exchanged to rearrange 
them so as to form a B meson moving at velocity v' . This rearrangement leads 
to a form-factor suppression, reflecting the fact that, as the velocities become 
more and more different, the probability for an elastic transition decreases. 
The important observation is that, in the limit mt —> c», the form factor can 
only depend on the Lorentz boost j = v ■ v' connecting the rest frames of the 
initial- and final-state mesons. Thus, in this limit a dimensionless probability 
function £,{v-v') describes the transition. It is called the Isgur- Wise functionll3. 
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In the HQET, which provides the appropriate framework for taking the limit 
nib —f cxD, the hadronic matrix element describing the scattering process can 
thus be written as 

— {B{v')\K,j^'b, \B{v)) = C(« • v') {y + v'T ■ (45) 

TUB 

Here b^ and 6^/ are the velocity-dependent heavy-quark fields of the HQET. 
It is important that the function ^(w • v') does not depend on rrif,. The factor 
1/tob on the left-hand side compensates for a trivial dependence on the heavy- 
meson mass caused by the relativistic normalization of meson states, which is 
conventionally taken to be 

{B{p')\B{p)) = 2mBv'^ (27r)3 6^{p- p') . (46) 

Note that there is no term proportional to {v ~ v')^ in (p5|). This can be seen 
by contracting the matrix element with {v — w')^, which must give zero since 

fbv = bv and bv'f — bv. 





t<tQ t = tQ t>tQ 

Figure 7: Elastic transition induced by an external heavy-quark current. 

It is more conventional to write the above matrix element in terms of an 
elastic form factor Fdiq'^) depending on the momentum transfer q^ = {p — p')'^: 

{B{v')\b^^b\B{v))^F,i{q')ip + pr, (47) 

where p^'' = niBV^''. Comparing this with (Eq), we find that 

i^ci(g') = e(« • «') , q^ ^ -2ml{v ■ v' - 1) . (48) 

Because of current conservation, the elastic form factor is normalized to unity 
at q^ — 0. This condition implies the normalization of the Isgur-Wise function 
at the kinematic point v ■ v' = I, i.e. for v = v': 

e(l) = l. (49) 

It is in accordance with the intuitive argument that the probability for an 
elastic transition is unity if there is no velocity change. Since for v = v' the 
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final-state meson is at rest in the rest frame of the initial meson, the point 
w • w' = 1 is referred to as the zero-recoil limit. 

The heavy-quark flavour symmetry can be used to replace the b quark in 
the final-state meson by a c quark, thereby turning the B meson into a D 
meson. Then the scattering process turns into a weak decay process. In the 
infinite-mass limit, the replacement by/ —^ Cyi is a symmetry transformation, 
under which the effective Lagrangian is invariant. Hence, the matrix element 

, ^ {D{v')\ c.,,j^by \B{v)) = ^v ■ v') {v + v'r (50) 

y/mBiTiD 

is still determined by the same function ^(w • v'). This is interesting, since in 
general the matrix element of a flavour-changing current between two pseu- 
doscalar mesons is described by two form factors: 

{D{v')\crb\B{v))^f+iq')ip+pr-f-{q'){p~pr- m 

Comparing the above two equations, we find that 

f±{q ) = ^{v-v), 

2y/mBmD 

q^ = m% + m\i — 2niBniD v ■ v' . (52) 

Thus, the heavy-quark flavour symmetry relates two a priori independent form 

factors to one and the same function. Moreover, the normalization of the 

Isgur-Wise function at w • w' = 1 now implies a non-trivial normalization of 

the form factors f±{q^) at the point of maximum momentum transfer, g^j^x = 

{■niB - tod)^: 

2 N mB±mD 
/±(9max = r, , • 53 

The heavy-quark spin symmetry leads to additional relations among weak 
decay form factors. It can be used to relate matrix elements involving vector 
mesons to those involving pseudoscalar mesons. A vector meson with longi- 
tudinal polarization is related to a pseudoscalar meson by a rotation of the 
heavy-quark spin. Hence, the spin-symmetry transformation c|,, -^ c*,|-telates 
B -^ D with B —f D* transitions. The result of this transformation istHi: 

, ^ {D*{v',e)\-Cy,rby \B{v)) = ze^'^^^ elv'^vp ^{v ■ v') , 

^ (D*{v',e)\cy>j^'j5by\B{v)) = \e*^ {vv' + I) - v'^' e* ■ v\c{v ■ v') , 



y/mBniD 



(54) 
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where e denotes the polarization vector of the D* meson. Once again, the 
matrix elements are completely described in terms of the Isgur-Wise function. 
Now this is even more remarkable, since in general four form factors, V{q'^) for 
the vector current, and Ai{q'^), i = 0,1,2, for the axial current, are required 
to pararnetrize these matrix elements. In the heavy-quark limit, they obey the 
relations O 



ms + Tod* 



av-v')^V{q')^Ao{q^)=A,{q^) 



{rriB + moY 



Ai{q^). 



2'm.B'niD-' V ■ v' . 



(55) 



Equations ( p2D and ( p5|) summarize the relations imposed by heavy-quark 
symmetry on the weak decay form factors describing the semileptonic decay 
processes B ^ D £D and B — > D*£ v. These relations are model-independent 
consequences of QCD in the limit where ?7T,f,, ttIc 2> Aqcd- They play a crucial 
role in the determination of the CKM matrix element |Vcb|. In terms of the 
recoil variable w = V^' ^ the differential semileptonic decay rates in the heavy- 
quark limit becomeLS: 



AV{B^Div) _ Gj 

dw 
dT{B^ D*iP) 



487r3 
G2 



|y,fc|2 {ms + ruD? ml {w^ ~ 1)^/2 ^^w) , 



dw; 



7^ l^chP {rriB - mo^ f ml, ^/w^~l {w + if 



Aw rn^ — 2wmBm,D' + mj-,. 



w + 1 



irriB - mo* 



e{w). (56) 



These expressions receive symmetry-breaking corrections, since the masses of 
the heavy quarks are not infinitely large. Perturbative corrections of order 
a"(7TiQ) can be calculated order by order in perturbation theory. A more 
difficult task is to control the non-perturbative power corrections of order 
(Aqcd/"^q)"- The HQET provides a systematic framework for analysing 
these corrections. For the case of weak-decag-jform factors, the analysis of 
the l/mq corrections was performed by LukeL^I. Later, Falk and the present 
author have analysed the structure of I/jtiq corrections for both meson and 
baryon weak decay form factorsE^I. We shall not discuss these rather technical 
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issues in detail, but only mention the most important result of Luke's analy- 
sis. It copcerns the zero-recoil limit, where an analoaae of the Ademollo-Gatto 
theorem im can be proved. This is Luke's theorem O, which states that the 
matrix elements describing the leading l/mg corrections to weak decay ampli- 
tudes vanish at zero recoil. This theorem is valid to all orders in perturbation 
theory IEItIS'LEI. Most importantly, it protects the B-^ D*£D decay rate from 
receiving first-order l/mg corrections at zero recoilllj. [A similar statement is 
not true for the decay B -^ D iv. The reason is simple but somewhat subtle. 
Luke's theorem protects only those form factors not multiplied by kinematic 
factors that vanish for v — v' . By angular momentum conservation, the two 
pseudoscalar mesons in the decay B —^ D £9 must be in a relative p wave, and 
hence the amplitude is proportional to the velocity \vo\ of the D meson in the 
B-meson rest frame. This leads to a factor (w^ — 1) in the decay cate. In such 
a situation, kinematically suppressed form factors can contribute Ell.] 



3.2 Short- Distance Corrections 

In section 0, we have discussed the first two steps in the construction of the 
HQET. Integrating out the small components in the heavy-quark fields, a non- 
local effective action was derived, which was then expanded in a series of local 
operators. The effective Lagrangian obtained that way correctly reproduces 
the long-distance physics of the full theory (see Fig. H). It does not contain the 
short-distance physics correctly, however. The reason is obvious: a heavy quark 
participates in strong interactions through its coupling to giuons. These gluons 
can be soft or hard, i.e. their virtual momenta can be small, of the order of the 
confinement scale, or large, of the order of the heavy-quark mass. But hard 
gluons can resolve the spin and flavour quantum numbers of a heavy quark. 
Their effects lead to a renormalization of the coefficients of the operators in 
the HQET. A new feature of such short-distance corrections is that through 
the running coupling constant they induce a logarithmic dependence on the 
heavy-quark massEj. Since as {"mQ ) is small, these effects can be calculated in 
perturbation theory. 

Consider, as an example, the matrix elements of the vector current V = 
qji^Q. —Ln QCD this current is partially conserved and needs no renormal- 
izationEj. Its matrix elements are free of ultraviolet divergences. Still, these 
matrix elements have a logarithmic dependence on toq from the exchange of 
hard gluons with virtual momenta of the order of the heavy-quark mass. If one 
goes over to the effective theory by taking the limit mq — > oo, these logarithms 
diverge. Consequeatly, the vector current in the effective theory does require 
a renormalizationtZI. Its matrix elements depend on an arbitrary renormaliza- 
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tion scale /i, which separates the regions of short- and long-distance physics. 
If n is chosen such that Aqcd ^ /i ^ mq , the effective coupling constant in 
the region between /i and uiq is small, and perturbation theory can be used 
to compute the short-distance corrections. These corrections have to be added 
to the matrix elements of the effective theory, which contain the long-distance 
physics below the scale fj,. Schematically, then, the relation between matrix 
elements in the full and in the effective theory is 

(y(mQ) )qcd = CoimQ,fi) (Vo(m)>hqet + — — — (V'i(m))hqet + .■• , 

(57) 
where we have indicated that matrix elements in the full theory depend on niQ , 
whereas matrix elements in the effective theory are mass- independent, but do 
depend on the renormalization scale. The Wilson coefRcients Ci{mQ, ^) are 
defined by this relation. Order by order in perturbation theory, they can 
be computed from a comparison of the matrix elements in the two theories. 
Since the effective theory is constructed to reproduce correctly the low-energy 
behaviour of the full theory, this "matching" procedure is independent of any 
long-distance physics, such as infrared singularities, non-perturbative effects, 
the nature of the external states used in the matrix elements, etc. 

The calculation of the coefficient functions in perturbation theory uses the 
powerful methods of the renormalization group. It is in principle straightfor- 
ward, yet in practice rather tedious. A comprehensive discussion of most of the 
existing calculations of short-distance corrections in the HQET can be found 
in Ref. 24. 



3.3 Model- Independent Determination of \Vcb\ 

We will now discuss some of the most important applications and tests of 
the above formalism in the context of semileptonic decays of B mesons. A 
model-independent determination of the CKM matrix element \Vcb\ based on 
heavy-quark symmetry can be obtained by J3aeasuring the recoil spectrum of 
D* mesons produced in 5 ^ D*iP decays LJ. In the heavy-quark limit, the 
differential decay rate for this process has been given in (bq) . In order to allow 
for corrections to that limit, we write 

dr(B^D*iD) G% . ,2 3 r-^ -, ,,2 



Aw m\ — 2w mBtnu* + m\). 



1 + 



1 {niB - niD'Y 
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\V,b?T\w), 



(58) 

where the hadronic form factor T{w) coincides with the Isgur-Wise function 
up to symmetry-breaking corrections of order as{'mQ) and Aqcd/'^q- The 
idea is to measure the product \Vcb\ ^{w) as a function of w, and to extract 
|ycfc| from an extrapolation of the data to the zero- recoil point w — \^ where 
the B and the D* mesons have a common rest frame. At this kinematic point, 
heavy-quark symmetry helps us to calculate the normalization T(\) with small 
and controlled theoretical errors. Since the range of w values accessible in this 
decay is rather small (1 < w < 1.5), the extrapolation can be done using an 
expansion around w = 1: 



T{w)^T[l) l-g^{w-l) + c{w-lf 



(59) 



The slope g^ and the curvature c are treated as parameters. 



0.05 




Figure 8: CLEO data for the pr(»4uct \Vci,\T{w), as extracted from the recoil spectrum in 
B — > D*iP decaysE3- The line shows a linear fit to the data. 

Measurements ofJjne recoil spectrum have been performed first by the 
ARGUS EJ and CLEOEj Collaborations in expemnents operating at the T(4al 
resonance, and more recently by the ALEPH eI, DELPHI and OPAL 
Collaborations at LEP. As an example. Fig. || shows the data reported by 
the CLEO Collaboration. The results obtained by the various experimental 
groups are summarized in Table m. In the first analyses, the curvature term in 
( p9| ) was omitted, and the data were fitted with a linear forin-£acjfcor. Later, 
the effect of a non-zero curvature has been taken into account Bej. It can be 
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shown in a model-independent way that the shape of thfyiprni factor is highly 
constrained by analyticity and unitarity requirementsEj^Ell. In particular, the 
curvature at w = 1 is strongly correlated with the slope of the form factor. For 
the values of ^^ given in Table |l|, one obtains a small positive curvature, so 
that the results obtained from a linear fit are very little affected by including 
a curvature term. The weighted average of the experimental results is 



|14fc|^(l) = (34.1 ± 1.4) X 10-3^ 
f = 0.80 ± 0.09 . 



(60) 



Table 1: Values for IVcbl-^Cl) (in units of 10 '') and q^ extracted from measurements of 
the recoil spectrum in B — > D*iu decays. 



Reference 


Method 


|y,6|.F(i)(io-3) 


-^9 

Q 


ARGU^ 


Linear Fit 


38.8 ±4.3 ±2.5 


1.17 ±0.22 ±0.06 


CLEOE^L 


Linear Fit 


35.1 ±1.9 ±2.0 


0.84 ±0.12 ±0.08 


alepheI 


Quadratic Fit 


32.0 ±2.1 ±2.0 


0.37 ±0.26 ±0.14 


DELPHI 


Linear Fit 


31.9 ±1.8 ±1.9 


0.31 ±0.17 ±0.08 


Linear Fit 


35.0 ±1.9 ±2.3 


0.81 ±0.16 ±0.10 


OPALy 


Quadratic Fit 


32.8 ±1.9 ±2.2 


0.55 ±0.24 ±0.05 




Linear Fit 


32.5 ±1.7 


0.42 ±0.17 



Heavy-quark symmetry implies that the general structure of the symmetry- 
breaking corrections to the form factor at zero recoil isO 



^(1)=77A 1 + Ox 



A 



QCD 



rriQ 



const X 



A2 



77^(l + ,5i/™2), (61) 



where tja is a short-distance correction arising from the finite renormalization 
of the flavour-changing axial current at zero recoil, and (Si/m^ parametrizes 
second-order (and higher) power corrections. The absence of [first-order power 
corrections at zero recoil is a consequence of Lukels-ttiepremlla. The one-loop 
expression for t^a has been known for a long time BUS: 

Gts [M) / rrib + mc m^ ^ 
TT \mi, — rric rric 3^ 

The scale M in the running coupling constant can be fixed by adopting the 
prescription of Brodsky, Lepage and Mackenzie (BLM) E3, where it is iden- 
tified with the average virtuality of the gluon in the one-loop diagrams that 



VA 



1 



0.96. 



(62) 
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contribute to tja- If asiA-I) is defined in the MS scheme, the resuh isL3 AI sa 
O.bly/m^jrih. Several estimates of higher-order corrections to r]A have been 
discussed. A renormalization-group resummation of logarithms of the 
(as ln?7ih/?7ic)", as{as\nmi,/mc)" and mc/mi,{ash\mb/mcy^ leads toEJ' 
rjA ~ 0.985. On the other hand, a resummation of "renormalon-chain" con- 
tributions of the form /3g ~^p^, where /3o = 11 — |n/ is the first coefficient of 



3' 

the QCD /3-function, gives l3 tia ~ 0.945. Using these partial resummations 
to estimate the uncertainty gives tja = 0.965 ± 0.020. Recently, Czarnecki has 
improved this estimate by calculating tja at two-loop orderEj. His result, 

77^ = 0.960 ±0.007, (63) 



is in excellent agreement with the BLM-improved one-loop expression 
Here the error is taken to be the size of the two-loop correction. 

The analysis of the power corrections is more difficult, since it cannot 
rely on perturbation theory. Three approaches have been discussed: in the 
"exclusive approach" , all l/rrig operators in the HQET are classified and their 
matrix elements estimated, leading toEJo 5i/m2 = — (3± 2)%; the "inclusive 
apprp^ch" has been used to derive the bound ^i/m^ < —3%, and to estimate 



thatL3^i/m2 — —(7 ±3)%; the "hybrid approach" combines the virtues of tt 



l/n 

former two to obtain a more restrictive lower bound on (5i /m^ • This leads to^ 

,5i/,„2 =-0.055 ±0.025. (64) 

Combining the above results, adding the theoretical errors linearly to be 
conservative, gives 

.F(l) = 0.91 ±0.03 (65) 

for the normalization of the hadronic form factor at zero recoil. Thus, the 
corrections to the heavy-quark limit amount to a moderate decrease of the 
form factor of about 10%. This can be used to extract from the experimental 
result (pOh the model-independent value 

\Vcb\ = (37.5 ± 1.5cxp ± 1.2th) X 10-3 . (66) 

3.4 Bounds and Predictions for q^ 

The slope parameter g^ in the expansion of the physical form factor in (59) 



differs from the slope parameter g^ of the Isgur-Wise function by corrections 
that violate the heavy-quark syjnmetry. The short-distance corrections have 
been calculated, with the result Ej 

^2 = e2^(0.16±0.02) + O(l/mQ). (67) 
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Bjorken has shown that the slope of the Isgur-Wise function is related to the 
form factors of transitions of a ground-state heavy meson into excited states, in 

,-f — 1+ ^T. 3 + 



which the light degrees of freedom carry quantum numbers j 



or 



by 



a sum rule which is an expression of quark-hadron duality: in the heavy-quark 
limit, the inclusive sum of the probabilities for decays into hadronic states is 
equal to the probability for the free quark transition. ^JL, qjos, normalizes the 
latter probability to unity, the sum rule takes the formE 



1 = 



1 



I^MI' 



I 



\e>iw)\^ 



+ (u;-l)|2^|r(™)Mp + (u; + l)2^|r(;)Mp|+0[(u;-l)2], 

(68) 

where l,m,n label the radial excitations of states with the same spin-parity 
quantum numbers. The terms in the first line on the right-hand side of the 
sum rule correspond to transitions into states containing light constituents with 
quantum numbers J^ = 5 • The ground state gives a contribution propor- 
tional to the Isgur-Wise function, and excited states contribute proportionally 
to analogous functions ^'^'-'(w). Because at zero recoil these states must be or- 
thogonal to the ground state, it follows that C^'-'(l) = 0, and the corresponding 
contributions to (pS]) are of order {w — 1)^. The contributions in the second 



line correspond to transitions into states with j^ = i or 



3 + 



Because of 



the change in parity, these are p-wave transitions. The amplitudes are propor- 
tional to the velocity \vf\ = {vJ^ — 1)'^/'^ of the final state in the rest frame of 
the initial state, which explains the suppression factor {w — 1) in the decay 
probabilities. The fimc-tions tj{w) are the analogues of the Isgur-Wise function 
for these transitions t21l. Transitions into excited states with quantum numbers 
other than the above proceed via higher partial waves and are suppressed by 
at least a factor [w — 1)^. 

For w = 1, eq. ( pq ) reduces to the normalization condition for the Isgur- 
Wise function. The Bjorken sum rule is obtained by expanding in powers of 
[w — 1) and keeping terms of first order. Taking iiitQ_a£C0unt the definition of 



the slope parameter, ^'(1) = - 



'^, one finds that 



-Q' 



'1/2 



100 



101 



(1)P + 2E 1-3/1(1)1' 



1 
>4 



(69) 



Notice that the lower bound is due to the prefactor i(w -I- 1) of the first 
term in (pq) and is of purely kinematic origin. In the analogous sum rule for 
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Aq baryons, this factor is absent, and consequently the slope parameter, 
the baryon Isgur-Wise function is only subject to the trivial constraint oE2: 
g^>0. 

Voloshin has derived another sum rule involving the form factors for tran 
sitions into excited states, which is the analogue of the jkiotical sum rule" for 
the dipole scattering of light in atomic physics. It reads 
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rrt n 

where Ej are the excitation energies relative to the mass uim of the ground- 
state heavy meson. This relation can be combined with the Bjorken sum rule 
to obtain an upper bound for the slope parameter g^: 

1 mM^mg 

where -Emin denotes the minimum excitation energy. In the quark model, one 
expects that i?min ~ f^M — "t-q, and one may use this as an estimate to obtain 
£>2 < 0.75.| 

The above discussion of the sum rules ignores renormalization effects. Both 
perturbative and non-perturbative corrections to (p9h and {uW can be incorpo- 
rated using the OPE, where one introduces a momentum scale /i ~ few x Aqcd 
large enough for as(/i) and power corrections of order (Aqcd//^)" to be small, 
but otherwise . a s..srnall as possible sa_as to suppress the contributions from 
excited states E£ai23. The result isE£Zl e^in(M) < Q^ < ftnaxC^); where the 
boundary values are shown in Fig. H as a function of the scale /i. Assuming 
that the OPE works down to values /i w 0.8 GeV, one obtains rather tight 
bounds for the slope parameters: 

0.5<£»2<o.8, 0.5<g2<i.i. (72) 

The allowed region for g^ has been increased in order to account for the un- 



known l/mq corrections in the relation (67). The experimental result given 



in (60) falls inside this region. 

These bounds compare well with theoretical predictions for the slope pa- 
rameters. QCD sum rules have been used to calculate the slope of the Isgur- 
Wise function. The results obtained by various authors are g^ = 0.84 ± 0.02 



'^Strictly speaking, the lowest excited state contributing to the sum rule is D + vr, which 
has an excitation-energy spectrum with a threshold at m^. However, this spectrum is broad, 
so that this contribution will not invalidate the upper bound for ^■^ derived here. 
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Figure 9: Bounds for the slope parameter (?• following from the Bjorken and Voloshin sum 

rules. 



), 0.7 ±0.1 (presei^uthorM), 0.70 ± 0.25 (Blok and Shif- 
). The UKQCD Collaboration has nm 



(Ba^iLet al 

manEB), and 1.00 ± 0.02 (Narison 

sented a lattice calculation of the slope of the form factor T{w) , yielding 

g^ = O-Q^qJIq^. We stress that the sum-rule bounds in JT^ ) are largely model 

independent; model calculations in strong disagreement with these bounds 

should be discarded. 



3. 5 Measurements of B ^ D* 
Heavy- Quark Symmetry 



id B ~f D £ ly Form Factors and Tests of 



We have discussed earlier in this section that heavy-quark symmetry implies 
relations between the semileptonic form factors of heavy mesons. They receive 
symmetry-breaking corrections, which can be estimated using the HQET. The 
extent to which these relations hold can be tested experimentally by comparing 
the different form factors describing the decays B -^ D^*^iD at the same value 
of w. 

When the lepton mass is neglected, the differential decay distributions 
in _B ^ D*ii? decays can be parametrized by three helicity amplitudes, or 
equivalently by three independent combinations of form factors. It has been 
suggested that a goocLcliaice for three such quantities should be inspired by 
the heavy-quark limit C3'lli3. One thus defines a form factor /i^i(w), which up 
to symmetry-breaking corrections coincides with the Isgur-Wise function, and 
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two form-factor ratios 

Riiw) 

R2{W) 



1- 



(ms +mD*y 



{rriB +mD*y 



A2{q^) 
Ai{q^) 



(73) 



The relation between w and q^ has been given in (Bq) . This definition is such 
that in the heavy-quark hmit Ri{w) — R2{w) = 1 independently of w. 

To extract the functions hAi{w)^ Ri{w) and R2{w) from experimental data 
is a complicated task. However, HQET-based calculations suggest that the w 
dependence of the fopiL-factor ratios, which is induced by symmetry-breaking 



effects, is rather mildEij. Moreover, the form factor hAiiw) is expected to have 



a nearly linear shape over the accessible w range. This motivates to introduce 
three parameters g\^, Ri and R2 by 



hAi{w) ^ T{1) [1 ~ g\^{w ~ 1) 



(74) 



where T{1) = 0.91 ± 0.03 from (||). The CLEO Collaboration has extracted 
these three ppiameters from an analysis of the angular distributions in B — > 
D*£p decaystiil. The results are 



ell =0.91 ±0.16, i?i = 1.18 ±0.32, i?2 =0.71 ±0.23. 



(75) 



Using the HQET, one obtains an essentially model-independent prediction for 
the symmetry-breaking corrections to i?i, whereas the corrections to R2 are 
somewhat model dependent. To good approximations 



Ri ^ 1 + 



R2 



1-K- 



37r 

A 

2mr 



A 
2mc 

0.8 ±0.2 



1.3 ±0.1, 



(76) 



with K w 1 from QCD sum rules l1^ . 

defined iiLXE8|). Theoretical calculations EiE 

analyses Ealjsuggest that A w 0.45-0.65 GeV is the appropriate value to be 

used in one-loop calculations. A quark-rnndel calculation of i?i and R2 gives 



e A is the "binding energy" as 
^^^ as well as phenomenological 



results similar to the HQET predictions till: i?i « 1.15 and R2 w 0.91. The 
experimental data confirm the theoretical prediction that i?i > 1 and ii'2 < 1, 
although the errors are still large. 
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There is a model-independent relation between the three parameters de- 
termined from the analysis of angular distributions aud the slope parameter 
p^ extracted from the semileptonic spectrum. It readsl£j 

Q\i~Q''^]:{Rl-l) + ,7r^^ ^(1-^2). (77) 

o{mB — mD*) 

The CLEO data give 0.07 ± 0.20 for the difference of the slope parameters on 
the left-hand side, and 0.22 ± 0.18 for the right-hand side. Both values are 
compatible within errors. 

More recently, heavy-quark symmetry has also been tested by comparing 
the form factor J-{w) in B -^ D*£ v decays with the corresponding-form factor 
Q{'w) governing B ^ D iD decays. The theoretical predictionEZl'ti3 

^ = 1.08 ±0.06 (78) 

compares well with the experimental results for this ratio: 0.99 ±0.19 reported 



by the CLEO CoUaborationE!!, and 0.87±0.30 reported by the ALEPH CoUab- 
orationE^. In these analyses, it has also been tested that within experimental 
errors the shape of the two form factors agrees over the entire range of w values. 
The results of the analyses described above are very encouraging. Within 
errors, the experiments confirm the HQET predictions, starting to test them 
at the level of symmetry-breaking corrections. 

3.6 Decays to Charmless Final States 

For completeness, we will discuss briefly semileptonic _B-meson decays into 
charmless flnal states, although heavy-quark symmetry does not help much 
in the analysis of these processes. Recently, the CLEO Collaboration has 
reported a first signal for the exclusive semileptonic decay modes B ^ Tr£D 
and B —f plD. The underlying quark process for these transitions is b ^ uiD. 
Thus, these decays provide information on the strength of the CKM matrix 
element Vub- The observed branching fractions areu: 

B{B -^ttIv) = (1.8 ± 0.5) X 10"" , 

B{B^ piv) = (2.5l°J) X lO"-* . (79) 

The theoretical description of these heavy-to-light {h — > u) decays is more 
model dependent than that for heavy-to-heavy [b -^ c) transitions, because 
heavy-quark symmetry does not help to constrain the relevant hadronic form 
factors. A variety of calculations for such form factors exists, based on QCD 
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Table 2: Values for IVuij/Vcbl extracted from the CLEO measurement of exclusive semilep- 

tonic B decays into charmless final states, taking IVcbl = 0.040. The first error quoted is 

experimental, the second (when available) is theoretical. 



Method 




B 



B ^ plv 



Sum Rules 



Nariaon ti 



0.177 ±0.025 ±0.001 

0.117±0.016± 0.012 

0.095 ±0.013 



0.064;°°J° ± 0.003 
0.090l[j;J]}6 ± 0.015 



Lattice QCD 



0.115 ±0.016j:o|]n 
0.094 ±0.013 ±0.023 



Pert. QCD 



0.060 ±0.008 



Quark Models 



ISGW2 
Melikhovt^ 



0.098 ±0.014 

0.098 ±0.014 

0.086 ±0.012 

0.099 ±0.014 ±0.014 



QAf^-, +0.010 

0.054l!J""f' 
0.083 



-0.007 

+0.013 

'0.015 



0.097tg|^i^ ± 0.014 



sum rules, lattice gauge theory, perturbative QCD, or quark models. Table || 
contains a summary of values extracted for the ratio \Vub/Vcb\ from a selection 
of such calculations. Some approaches are more consistent than others in 
that the extracted values are compatible for the two decay modes. With few 
exceptions, the results lie in the range 



Vub 



Vcb 



0.06-0.11 . 



(80) 



cxcl 



which is in good agreement with the measurement of |Vu&| obtained froiiLthe 
endpoint region of the lepton spectrum in inclusive semileptonic decaysuu: 



Kb 

Vcb 



0.08 ±0.01cxp± 0.02th 



(81) 



incl 



Clearly, this is only the first step towards a more reliable determination of 
iKtbl; yet, with the discovery of exclusive b —^ u transitions an important mile- 
stone has been met. Efforts must now concentrate on more reliable methods 
to determine the form factors for heavy-to-light transitions. Some new ideas in 

this dicfiction have been discussed recently. They are based on lattice calcula- 

-sdir- 



tionst23, analyticity constrain^ 
for heavy-to-heavy transitions 



13C 



132 



131 



or variants of the form-factor relations 
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4 Inclusive Decay Rates and Lifetimes 

Inclusive decay rates determine the probability of the decay of a particle into 
the sum of all possible final states with a given set of global quantum numbers. 
An example is provided by the inclusive scmilcptonic decay rate of the B 
meson, T{B —^ X iD), where the final state consists of a lepton-neutrino pair 
accompanied by any number of hadrons. Here we shall discuss the th^Qi:^t,ic£^l 
description of inclusive decays of hadrons containing a heavy quark 
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From the theoretical point of view, such decays have two advantages: first. 



bound-state effects related to the 



la. 



the heavy quark inside the hadronElHIElil) can be accounted for in a systematic 



state (such as the "Fermi motion" of 



way using the heavy-quark expansion; secondly, the fact that the final state 
consists of a sum over many hadronic channels eliminates bound-state effects 
related to the properties of individual hadrons. This second feature is based on 
the hypothesis of quark-hadron duality, which is an important concept in QCD 
phenomenology. The assumption of duality is that cross sections and decay 
rates, which are defined in the physical region (i.e. the region of time-like 
momenta), are caJmlable in QCD after a "smearing" or "averaging" procedure 



has been appliedE13. In semileptonic decays, it is the integration over the lepton 
and neutrino phase space that provides a smearing over the invariant hadronic 
mass of the final state (so-called global duality). For non-leptonic decays, on 
the other hand, the total hadronic mass is fixed, and it is only the fact that 
one sums over many hadronic states that provides an averaging (so-called local 
duality) . Clearly, local duality is a stronger assumption than global duality. It 
is important to stress that quark-hadron duality cannot yet be derived from 
first principles; still, it is a necessary assumption for many applications of 
QCD. The validity of |global duality has been tested experimentally using data 
on hadronic r decays till. Some more formal attempts to address the problem 
of quark-hadron duality can be found in Refs. 145,146. 

Using the optical theorem, the inclusive decay width of a hadron Hb con- 
taining a b quark can be written in the form 

r{Hb^X) = ^lm{Hb\T\Hb), (82) 

where the transition operator T is given by 

T = ifd^xT{C,six),CeffiO)}. (83) 

Inserting a complete set of states inside the time-ordered product, we recover 
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the standard expression 



2to//, ^ 



for the decay rate. For the case of semileptonic and non-leptonic decays, Ccs 
is the effective weak T agranrriayi giypn in (m, which in practice is corrected 



for short-distance effects c3'E3'tl3iilj arising from the exchange of gluons with 



virtuahties between mw and nib- If some quantum numbers of the final states 
X are specified, the sum over intermediate states is restricted appropriately. In 
the case of the inclusive semileptonic decay rate, for instance, the sum would 
include only those states X containing a lepton-neutrino pair. 

— » — n — » n — » — — ► — » — • — • — 







bgsa^,j,Gl-"'b 



Figure 10: Perturbative contributions to the transition operator T (left), and the corre- 
sponding operators in the OPE (right). The open squares represent a four-fermion interaction 
of the effective Lagrangian £effi while the black circles represent local operators in the OPE. 

In perturbation theory, some contributions to the transition operator are 
given by the two-loop diagrams shown on the left-hand side in Fig. ^ Because 
of the large mass of the b quark, the momenta flowing through the internal 
propagator lines are large. It is thus possible to construct an OPE for the 
transition operator, in which T is represented as a series of local operators 
containing the heavy-quark fields. The operator with the lowest dimension, 
d = 3, is 56. It arises by contracting the internal lines of the first diagram. 
The only gauge-invariant operator with dimension 4 is bipb; however, the 
equations of motion imply that between physical states this operator can be 
replaced by m^bb. The first operator that is different from bb has dimension 5 
and contains the gluon field. It is given by bggO'^i.iyG^'^b. This operator arises 
from diagrams in which a gluon is emitted from one of the internal lines, such 



as the second diagram shown in Fig. 10. For dimensional reasons, the matrix 
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elements of such higher-dimensional operators are suppressed by inverse powers 
of the heavy-auaiJi_mass. Thus, any inclusive decay rate of a hadron Hh can 
be written aslHJT^^ 

^^^^^^f^ = ^^^Y^^^^^" + 'i ^^2 + ---|^ (85) 

where the prcfactor arises naturally from the loop integrations, c^ are calcu- 
lable coefficient functions (which also contain the relevant CKM matrix ele- 
ments) depending on the quantum numbers / of the final state, and {0)h are 
the (normalized) forward matrix elements of local operators, for which we use 
the short-hand notation 

{0)h = y^ {Ht\ O \H,) . (86) 

In the next step, these matrix elements are systematically ex panded in 
powers of l/rn-b, using the technology of the HQET. The result is EJ'P4P^ 

,2/-fr \ _ ,,2 



{hh)n^l-^^^^^^^^^P^^ + 0{l/ml)., 



2ml 

{bgs'yt..G'"'b)H - 2fi%{Hb) + 0(l/mb) , (87) 

where we have defined the HQET matrix elements 

^il{H,) = -^ {Ht{v)\ K (iDf K \Hb{v)) , 
2mHt, 

f^UHb) = TT^ {Hb{v)\ k^<J^.G^''b, \H,{v)) . (88) 

2mHt 2 

Here {iD)^ = {iv ■ D)^ — {iD)^; in the rest frame, this is the square of the 
operator for the spatial momentum of the heavy quark. Inserting these results 
into (H) yields 

(89) 
It is instructive to understand the appearance of the "kinetic energy" contri- 
bution /xj, which is the gauge-covariant extension of the square of the 6-quark 
momentum inside the heavy hadron. This contribution is the field-theory ana- 
logue of the Lorentz factor (1 — v^Y^^ ~ 1 — fc^/2m^, in accordance with the 
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fact that the hfetime, r — l/F, for a moving particle increases due to time 
dilation. 

The main result of the heavy-quark expansion for inclusive decay rates 
is the observation that the free quark decay (i.e_the parton model) provides 
the first term in a systematic 1/mt, expansion t^. For dimensional reasons, 
the corresponding rate is proportional to the fifth power of the 5-quark mass. 
The non-perturbative corrections, which arise from bound-state effects inside 
the B meson, are suppressed by at least two powers of the heavy-quark mass, 
i.e. they are of relative order (AqcD/n^b)^- Note that the absence of first- 
order power corrections is a consequence of the equations of motion, as there 
is no independent gauge-invariant operator of dimension 4 that could appear 
in the OPE. The fact that bound-state effects in inclusive decays are strongly 
suppressed exj:)lai,ps_fi posteriori the success of the parton model in describing 
such processes 
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The hadronic matrix elements appearing in the heavy-quark expansion 
( p9| ) can be determined to some extent from the known masses of heavy hadron 
states. For the B meson, one finds that 

Ai^(S) = -Ai = (0.3 ± 0.2) GeV^ , 

fil.{B) = 3A2 w 0.36 GeV^ , (90) 

where Ai and A2 are the parameters appearing in the mass formula (p2). For 
the ground-state baryon A;,, in which the light constituents have total spin 
zero, it follows that 

MG(Afc) = 0, (91) 

while the matrix element /i^(Ab) obeys the relation 

(mA, -mAj -(ms -mi)) = nl{B) - nl{Kb) ( ^;^ " ^^^ j + '^(1/"^q) ' 

(92) 
where mB and ruD denote the spin-averaged masses introduced in connection 
with (E3). With the value of toa^ given in (EQ), this leads to 

^il{B) - ^il{Kb) = (0.01 ± 0.03) GeV^ . (93) 

What remains to be calculated, then, is the coefhcient functions c^ for a given 
inclusive decay channel. We shall now discuss some of the most important 
applications of this general formalism. 
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4-.1 Determination of \Vcb\ from Inclusive Semileptonic Decays 

The extraction of | Vcb \ from the inclusive semileptonic decay rate of B meaansJi 



based on the general expression (p9) , with the short-distance coefficientsE£i 
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„SL _ I-,/ |2 



Vcl 



1 - 8x^ + Sx" - a;"* - 12a;* In a;^ + 0{a, 



„SL 



-6\V,b\\l-xY- (94) 



Here x = mc/mi,, and mt and rric are the pole^nasses of the b and c quarks, 
defined to a giveii order in perturbation theoryc3. The 0{as) terms in cf" are 
knowp-^pouptly t£3, while only partial calculations of higher-order corrections 
The theoretical uncertainties in this determination of \Vcb\ are 



IM for the central 
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quite different from those entering the analysis of exclusive decays. The main 

sources are the dependence on the heavy-quark masses, unknown higher-order 

perturbative corrections, and the assumption of global quark-hadron duality. 

A consetifative estimate of the total theoretical error on the extracted value of 

\Vcb\ isN S\Vcb\/\Vcb\ « 10%. Taking the result of Ball et al 

value, and using tb = (1.60 ±0.03) ps for the average i3-meson lifetime 

find 

|yefc| = (0.040 ± 0.004) f^^j = (40 ± lexp ± 4th) X 10-3 . (95) 

In the last step, we have used BgL = (10.8±0.5)% for the semileptonic branch- 
ing ratio of B mesons (see below). The value of \Vcb\ extracted from the in- 



clusive semileptonic width is in excellent agreement with the value in (66) 
obtained from the analysis of the exclusive decay B -^ D*£D. This agreement 
is gratifying given the differences of the methods used, and it provides an in- 
direct test of global quark-hadron duality. Combining the two measurements 
gives the final result 

I Vcb I =0.039 ±0.002. (96) 

After Vud and Ku, this is the third-best known entry in the CKM matrix. 

4-2 Semileptonic Branching Ratio for Decays into r Leptons 

Semileptonic decays of B mesons into r leptons are of particular inmartance, 
since they are sensitive probes of physics beyond the Standard ModelE£3. From 
the theoretical point of view, the ratio of the semileptonic rates (or branching 
ratios) into r leptons and electrons can be calculated reliably. This ratio is 
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independent of the factor ml, the )aaijjixi|[^ic_paranieter Ai, and CKM matrix 
elements. To order 1/m^, one finds 



156,159 



-^ —^ = fix,, xr) + ^ gixc, Xr) = 0.22 ± 0.02 , 97 

where / and g are calculable coefficient functions depending on the mass ratios 
Xc = rricl'm^i and x-^ = mT/mf,, as well as on as{m},). Two new measurements 
of the semileptonic branching ratio of b quarks, B(& -^ Xjip^, have been re- 



ported by the ALEPH and OPAL Collaborations at LEPE2!J'm. The weightec 
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average is (2.68 ±0.28)%. Normalizing this result to the LEP average value 
B(5 -> X e Ve) = (10.95 ± 0.32)%, we obtain 

11^14^ = 0.245 ±0.027, (98) 

in good agreement with the theoretical prediction ( |97| ) for B mesons. 

4-. 3 Semileptonic Branching Ratio and Charm Counting 
The semileptonic branching ratio of B mesons is defined as 

E^r(B->x£z>) + rhad + r,arc' ^ ^ 

where Fhad and Praro are the inclusive rates for hadronic and rare decays, 
respectively. The main difficulty in calculating BgL is not in the semileptonic 
width, but in the non-leptonic one. As mentioned previously, the calculation 
of non-leptonic decay rates in the heavy-quark expansion relies on the strong 
assumption of local quark-hadron duality. 

Measurements of the semileptonic branching ratio have been performed 
by various experimental groups, using both model-dependent and model-inde- 
pendent analyses. The status of the results is controversial, as there is a 
discrepancy between low-energy measurements performed at the T(4s) reso- 
nance and high-energy measurements performed-al the Z^ resonance. The 



situation has been reviewed recently by Richmant^, whose numbers we shall 
use. The average value at low energies is Bsl = (10.23 ± 0.39)%, whereas 
high-energy measurements give Bsl(&) = (10.95 ± 0.32)%. The label (6) indi- 
cates that this value refers not to the B meson, but to a mixture of b hadrons 
(approximately 40% B", 40% B° , 12% B^, and 8% A^). Assuming that the 
corresponding semileptonic width Psl(^) is close to that of B mesons,f| we can 

''Theoretically, this is expected to be a very good approximation. 
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correct for this fact and find Bsl = (tb/t;,) BslI^) = (11.23 ± 0.34)%, where 
Tf, = (1.56 ±0.03) ps is the average hfetime corresponding to the above mixture 
of b hadrons. The discrepancy between the low- and high-energy measurements 
of the semileptonic branching ratio is therefore larger than three standard de- 
viations. If we take the average and inflate the error to account for this fact, 
we obtain 

Bsl = (10.80 ±0.51)%. (100) 

An important aspect in interpreting this result is charm counting, i.e. the 
measurement of the average number ric of charm hadrons produced per B 
decay. Theoretically, this quantity is given by 

nc = l + B{B^ Xcc.') - B{B -^ X„oc) , (101) 

where B(B — > Xcss') is the branching ratio for decays into final states con- 
taining two charm quarks, and B{B =i_.X_c) ~ 0.02 is the Standard Model 
branching ratip_£pr charmless decays ll£3~li£3. The average value obtained_a±, 



low energies isEr^ n^ — 1.12 ± 0.05, whereas high-energy measurements give 
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ric = 1.23 ± 0.07. The weighted average is 

nc = 1.16 ±0.04. (102) 

The naive parton model predicts that Bsl ~ 15% and ric ~ 1.2; however, it 
has been known for some time that perturbative corrections could change these 



results significantlyE£3. With the establishment of the heavy-quark expansion, 
the non-perturbative corrections to the parton model could be computed, and 
their effect turned out to be very small. This led Bigi |et-,al. to conclude that 



values Bsl < 12.5% cannot be accommodated by theoryE£3. Later, Bagan et al 



have completed the calculation of the 0{as) corrections including the effects-qf 
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the charm-quark mass, finding that they lower the value of BgL significantly 
Their original analysis has recently been corrected in an erratum. Here we 
shall present the results of an independent numerical analysis using the same 



theoretical input ll£3. The semileptonic branching ratio and ric depend on 
the quark-mass ratio mc/mh and on the ratio /i/TTif,, where /i is the scale 
used to renormalize the coupling constant as(/i) and the Wilson cocfhcients 
appearing in the non-leptonic decay rate. The freedom in choosing the scale 
/x reflects our ignorance of higher-order corrections, which are neglected when 
the perturbative expansion is truncated at order as- We allow the pole masses 
of the heavy quarks to vary in the range ttij, = (4.8 ± 0.2) GeV and mi, — rUc = 
(3.39 ± 0.06) GeV, corresponding to 0.25 < mJjni, < 0.33. The value of 
the difference {mi, — mc) is in accordance with m%. Non-perturbative effects 
appearing at order 1/to^ in the heavy-quark expansion are described by the 
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single parameter A2 in (p7|); the dependence on the parameter Ai is the same 
for all inclusive decay rates and cancels out in the preilicjtions for BgL and Uc 
For the two choices /i = rrn, and /j, = m(,/2, we obtainE 
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BgL — 



12.0 ±1.0%; At = "^b, 

10.9 ±1.0%; /x = m6/2, 

1.20 T 0.06; fi = mt, 

1.21 =F 0.06; n^rrn,/2. 



(103) 



The uncertainties in the two quantities, which result from the variation of 
mc/Viib in the range given above, are anticorrelated. Notice that the semilep- 
tonic branching ratio has a stronger scale dependence than jic- By choosing a 
low renormalization scale, values Bsl < 12.5% can easily be accompiodated. 
This is indeed not unnatural. Using the BLM scale-setting method EJ, it Jiap 
been estimated that /i > 0.32r7T.(, is an appropriate scale to use in this case 
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Figure 11: Theoretical prediction for the semileptonic branching ratio and charm counting 

as a function of the quark-mass ratio mc/mi, and the renormalization scale fi. The data 

points show the average experimental values obtained in low-energy (LE) and high-energy 

(HE) measurements, as discussed in the text. 

The combined theoretical predictions for the semileptonic branching ra- 
tio and charm counting are shown in Fig. O. They are compared with the 
experimental results obtained from low- and high-energy measurements. It 
has been argued that the combination of a low semileptonic branching ratio 
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and a law value of ric would constitute a potential problem for the Standard 



Modelll£j. However, with the new experimental and theoretical numbers, only 



for the low-energy measurements a discrepancy remains between theory and 



experiment. Note that, with (101), our results for Uc can be used to calculate 
the branching ratio B(i? — > Xccs')i which is accessible to a direct experimental 
determination. Our prediction of (22 ±6)% for this branching ratio agrees waL 
with the preliminary result reported by the CLEO Collaboration, which is 
B(B ^ X,E,0 = (23.9 ± 3. 
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Previous attempts to resolve the "problem of the semileptonic branching 
ratio" have focused on four possibilities: 

• It has been argued that the experimental value of ric may depend on 
model assumptions abnnt thp production of charm hadrons, which are 



sometimes questionable E£jEZj. 



• It has been pointed out that the asswiuptioD_Qf local quark- hadron duality 
could fail in non-leptonic B decayslll3'tIiHll3. If so, this will most likely 
happen in the channel b —f ccs, where the energy release, E = vriB — 
"T-x(ccs): is less than about 1.5 GeV. However, if one assumes that sizeable 
duality violations occur only in this channel, il-ia, impossible to improve 
the agreement between theory and experiment!] 
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• Another possibility is that higher-order corrections in the l/rrii, expan- 
sion, which were previously thought to be negligible, give a sizeable con- 
tribution. It has been argued that they could lower the semileptonic 
branching ratio Jajc, up to 1%, depending on the size of some hadronic 



matrix elementsllS^ 
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1. 



Finally, there is also the possibility to invoke New Physics lil3iiZ3. One 
may, for instance, consider extensions of the Standard Model with en- 
hanced flavour-changing neutral currents, such as 6 — > sg. The effect 
of such a contribution would be that both Bsl and Uc are reduced by a 
factor (1 -f 77B|^)-\ where 77 = (r,aro - r^are)/rsL. To obtain a sizeable 
decrease requires values 77 > 0.5, which are large (in the Standard Model, 
TrEtra/rsL ~ 0.2), but not excluded by current experiments. 



4.4 Lifetime Ratios of b Hadrons 

The heavy-quark expansion shows that the lifetimes of all hadrons containing 
a b quark agree up to non-perturbative corrections suppressed by at least two 
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powers of l/rrib. In particular, it predicts that 
r{B-) 



1 + Oil/ml), 

{1.00 ±0.01) + Oil/ml), 

r(Ab) f4iAb) ~ f^lJB) tiUB) 

Vm~ ^, ca^^ + Oil/m,) 

« 0.98 + O(l/mg), (104) 



t(BO) 

rjBs) 
r{Bd) 



where E£3 cq w 1.1, and we have used (90) and (P3|). The uncertainty in the 



value of the ratio T{Bs)/T{Bd) arises from unknown SU(3)-violating effects in 
the matrix elements of Bg mesons. The above theoretical predictions may be 
conmared with the average experimental values for the lifetime ratios, which 
areEi: 



^{B- 



rjBs) 
r{Bd) 

r(Afc) 
r(BO) 



= 1.06 ±0.04, 
= 0.98 + 0.07, 
= 0.78 ±0.04. (105) 



Whereas the lifetime ratios of the different B mesons are in good agreement 
with the theoretical prediction, the low value of the lifetime of the Af, baryon 
is surprising. 

To understand the atiiicture of lifetime differences requires to go further 
in the l/m;, expansion llZ3. Although at first sight it appears that higher- 
order corrections could be safely neglected given the smallness of the 1/m^ 
corrections, this impression is erroneous for two reasons: first, at order l/m^ 
in the heavy-quark expansion for non-leptonic decay rates there appear four- 
quark operators, whose matrix elements explicitly depend on the flavour of 
the spectator quark(s) in the hadron iJf,, and hence are responsible for lifetime 
differences between hadrons with different light constituents; secondly, these 
spectator effects receive a phase-space e nha ncement factor of 0(167r^) with 
respect to the leading terms in the 0PEE£3. This can be seen from Fig. [l2| , 
which shows that the corresponding contributions to the transition operator 
T arise from one- loop rather than two-loop diagrams. The presence of this 
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Figure 12: Spectator contributions to the transition operator T (left), and the correspond- 
ing operators in the OPE (right). Here F; denotes some combination of Dirac and colour 

matrices. 



phase-space enhancement factor leads to a pecuhar structure of the heavy- 
quark expansion for non-lcptonic rates, which may be displayed as follows: 



Toil + X2 



A, 



QCD 

mb 



/Aqcd 
X3[ —^ 

V mb 



167r^ 



/Aqcd 
2/3 -^ — 
V mb 



2/4 ^^^ — 
mb J 



(106) 



Here Fq is the free quark decay rate, and Xn and j/„ are coefficients of order 
unity. Since it is conceivable that the terms of order IGtt^ {-^qcd / m,b)^ could 
be larger than the ones of order (Aqcd /'^ib)^, it is important to include them 
in the predictions for non-lcptonic decay rates. Moreover, there is a challenge 
to calculate the hadronic matrix elements of the corresponding four-quark op- 
erators with high accuracy. 

In total, a set of four four-quark operators is induced by spectator effects. 
They are 



rpq 
rpq 

-'s-p 



&7p(l- 75)997^(1 -75)&, 

^(1-75)99(1 + 75)^, 

&7p(l - l5)taqqi''{l - l5)tab, 

6(1 -75)^099(1 +l5)tab, 



(107) 



where q is a light quark, and ta are the generators of colour SU(3). In most 
previous analyses of spectator effects the hadronic matrix elements of these op- 
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erators have been estimated making simplifying assumptions EZjni£3. For the 
matpx elements between _B-meson states the vacuum saturation approxima- 
tionCj has been assumed, i.e. the matrix elements of the four-quark operators 
have been evaluated by inserting the vacuum inside the current products. This 
leads to 

(^,1 0^_A \B,) = {B,\ Ol_p \B,) = flm\ , 

(^,1 T^_A \B,) = {B,\ Tl_p \B,) = , (108) 

where fs is the decay constant of the Bq meson, defined as 

(0 I q-ff^js b \Bqiv)) = ifBrnsy . (109) 

The vacuum saturation approximation has been criticized by ChernyakE£3, 
who estimates that the corrections to it can be as large as 50%. 

A model-independent analysis of spectator effects, which avoids assump- 
tions about hadronic matrix elements, can be performed if instead of (|108D one 



defines Il68 



{Bq\Ol^p\Bq) =B2flm\, 
{Bq\T^_j,\Bq) -£i/I,m|^, 
{Bq\Tl_^\Bq)=e2flml,- (HO) 

The values of the dimensionless hadronic parameters Bi and Si are currently 
not known; ultimately, they may be calculated using some field-theoretic ap- 
proach such as lattice gauge theory or QCD sum rules. The vacuum saturation 
approximation corresponds to setting Bi = 1 and e^ = (at some scale fi, where 
the approximation is believed to be valid). For real QCD, it is known that 

B, -0(1), £, = 0(1/7V,), (111) 

where Nc is the number of colours. Below, we shall treat Bi and Si (renormal- 
ized at the scale m;,) as unknown parameters. Similarly, the relevant hadronic 
matrix elements of the four-quark operators between A(,-baryon states can be 
parametrized by two parameters, B and r, where _B = 1 in the valence-quark 
approximation, in which the colour of the quark fields in the operators is iden- 
tified with the colour of the quarks inside the baryon. 
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Lifetime ratio for B and B^ 



The lifetimes of the charged and neutral B mesons differ because of two types of 
spectator effects iliuatrnted in Fig. O. They are referred to as Pauli interference 
and W exchangeE^"M3. In the operator language, these effects are represented 



by the hadronic matrix elements of the local four-quark operators given in 



(110). In fact, the diagrams in Fig. [13 can be obtained from those in Fig. 12 



by cutting the internal lines, which corresponds to taking the imaginary part 
in 




Figure 13: Pauli interference and W exchange contributions to the lifetimes of the B~ and 

the B'^ mesons. The spectator effect in the first diagram arises from the interference due to 

the presence of two identical u quarks in the final state. 



The explicit calculation of these spectator effects leads to 
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Arspcc(-Bg) — 



G 



1927r3 



Vcbr 16^ 



2 
2 JB 



JbI^b 



m-i 



(112) 



where 



Cb- 



-0.4Bi + 6.6ei , 



Cso 



-2.2ei + 2.4e2 



(113) 



Note the factor of 167r^ in ( 112 ), which arises from the phase-space enhance- 
ment of spectator effects. Given that the parton-model result for the total 
decay width is 



rtot(B)«3.7x 



llml 



\Vcb\\ 



1927r3 
we see that the characteristic scale of the spectator contributions is 

2 



47r 



2 Ib^t-b _ /27r/, 



V "16 



5%. 



(114) 



(115) 



The precise value of the lifetime ratio depends crucially on the size of the 
hadronic matrix elements. Taking Jb = 200 MeV for the decay constant of the 
B meson (see Ref. 24 and references therein), i.e. absorbing the uncertainty in 
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this parameter into the definition of Bi and Si, we find 
r{B-) 



Wi 



t{BO) 



1 + 0.03Bi - 0.70ei + 0.20£2 

1 + 0.05i?i - 0.75ei + 0.20£2 ■ (116) 



Here Bi and Si refer to a renornialization scale /i = mi,, whereas Bi and e^ 
refer to a low renornialization point /i — /ihad chosen such that as(/ihad) = 0.5. 
The most striking feature of this result is that the coefhcients of the colour- 
octet operators Ty-A and Tg-p are an order of magnitude larger than those 
of the colour-singlet operator Oy-A- As a consequence, the vacuum insertion 
approximation, which was adopted in Ref. 179 to predict that t{B^)/t{B'^) 
is larger than unity by an amount of order 5%, should not be trusted (not 
even at a low renornialization point such as /Xhad)- With Si of order 1/iVc, 
it is conceivable that the non-factorizable contributions dominate the result. 
Thus, without a detailed calculation of the parameters Si no reliable prediction 
can be obtained. Given our present ignorance about the true values of the 
hadronic matrix elements, we must conclude that even the sign of the sum of 
the spectator contributions cannot be predicted. A lifetime ratio in the range 
0.8 < t{B~)/t{B'^) < 1.2 could be easily accommodated by theory. 

In view of these considerations, the experimental fact that the lifetime 
ratio turns out to be close to unity is somewhat of a surprise. It implies a 
constraint on a certain combination of the colour-octet matrix elements, which 
reads ei — 0.3e2 = few %. 



Lifetime ratio for Bg and Bd 

The lifetimes of the two neutral mesons Bg and Bd differ because spectator 
effects depend on the flavour of the light quark, and moreover because the 
hadronic matrix elements in the two cases differ by SU(3) symmetry- breaking 
corrections. It is difficult to predict the sign ot.thij; net effect, but the magnitude 



cannot be larger than one or two per centllSaliZj. Hence 

riBs) 



r{Bd) 



= 1±0(1%), (117) 



denotes the average lifetime of the two i?s_si.ates, whose individual lifetimes 



which is consistent with the experimental value in (105). Note that t{Bs) 
denotes the average lifetime of the two Bg. 
are expected to differ by a sizeable amountt 

46 



174[l84 



Lifetime ratio for Ah and S" 



Although, as shown m (104), the hfctimc differences between heavy mesons and 
baryons start at order l/m,^ in the heavy-quarlc expansion, the main effects 
are expected to appear at order 1/m^. Therefore, for an estimate of the ratio 
t(A(,)/t(_B'') one needs the matrix elements of four-quark operators between 
baryon states. Very little is known about such matrix elements. Bjigi, et al. 
have adopted a simple non-relativistic quark model to conclude that 

r(Ab) 
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t{BO) 



0.90-0.95. (118) 



An even smaller lifetime difference has haen obtained by Rosnert£: 



A model-independent analysis gives 
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rm 



0.98 - 0.17£i + 0.2062 - (0.012 -I- 0.021B)r , (119) 



where B and r are expected to be positive and of order unity. Given the 
structure of this result, it seems difficult to explain the experimental value 
T(Ab)/r(_B") = 0.78 ± 0.04 without violating the constraint ei w 0.3e2 men- 
tioned above."! Essentially the only possibility is to have r of order 2-4 or so, 
as there are good theoretical arguments why B cannot be much larger than 
unity. On the other hand, in a constituent quark picture, r is the ratio of the 
wave functions determining the probability to find a light quark at the location 
of the b quark inside the At, baryon and the B meson, i.e. 

.J<^. (120, 

l»'5(tl)P 

and it is hard to see how this ratio could be much different from unity. 

In view of the above discussion, the observation of the short A^ lifetime 
remains a puzzle, whose explanation may lie beyond the heavy-quark expan- 
sion. If the current experimental value persists, one may have to question 
the validity of local quark-hadron duality, which is assumed in-t^e theoretical 
calculation of lifetimes and non-leptonic inclusive decay rates 
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5 Concluding Remarks 

We have presented a review of the theory and phenomenology of heavy-quark 
symmetry and its applications to the spectroscopy, weak decays and lifetimes 

■^Another constraint arises if one does not want to spoil the theoretical prediction for the 
semileptonic branching ratio. 
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of hadrons containing a heavy quark. The theoretical tools that allow us to 
perform quantitative calculations are the heavy-quark effective theory and the 
1/toq expansion. Our hope is to have convinced the reader that heavy-flavour 
physics is a rich and diverse area of research, which is at present characterized 
by a fruitful interplay between theory and experiments. This has led to many 
significant discoveries and developments. Heavy-quark physics has the poten- 
tial to determine many important parameters of the electroweak theory and 
to test the Standard Model at low energies. At the same time, it provides an 
ideal laboratory to study the nature of non-perturbative phenomena in QCD. 
The prospects for further significant developments in the field of heavy- 
flavour physics look rather promising. With the approval of the first asymmet- 
ric B factories at SLAG and KEK, with ongoing i3-physics programs at the 
existing facilities at Cornell, Fermilab and CERN, and with plans for future B 
physics at HERA-B and the LHC-B, there are i?cautiful times ahead of us! 
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